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We derive an effective hamiltonian in the Lowest Landau Level (LLL) that incorporates the effects 
of Landau-level mixing to all higher Landau levels to leading order in the ratio of interaction energy 
to the cyclotron energy, k = (e 2 /el)/uj c . We then transcribe the hamiltonian to the composite 
fcrmion basis using our hamiltonian approach and compute the effect of LL mixing on transport 
gaps. 
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I. INTRODUCTION 

The fractional quantum Hall effects (FQHE) 0,0] have 
demonstrated the existence of new states of matter for 
two-dimensional electron gases (2DEG's) in high mag- 
netic fields. When the fields are high enough that only a 
fraction of a Landau level (LL) is occupied, the system 
reorganizes itself at special fractions into incompressible 
ground states j|] with fractionally charged excitations. 
The current understanding of the FQHE is schematically 
expressed by saying that the quasiparticles are electrons 
dressed by an even number of units of statistical flux, 
called Composite Fermions (CF's) The FQHE of 

electrons is then understood as the integer quantum Hall 
effect of CF's. 

A fruitful approximation in the study of the FQHE 
has been the restriction of the single particle electronic 
wavefunctions to the LLL. This approximation becomes 
exact in the limit of very high magnetic fields (k = 
(e 2 /el)/uj c — > 0) and captures the essential physics under- 
lying the FQHE. However, when one wishes to make de- 
tailed comparisons between the theory and experiment, 
say in the study of transport gaps, one must seek cor- 
rections due to the nonzero value of k, that is, include 
the effects of Landau level mixing. Hitherto these effects 
have been studied numerically with the latest re- 

sults coming from a fixed-phase diffusion Monte Carlo 
study §. 

In this paper we present an analytical calculation of 
the effect of LL-mixing in two steps. First we deduce an 
effective hamiltonian within the LLL which incorporates 
the leading contribution in k due to higher LL's. This 
hamiltonian is written in terms of electronic coordinates. 
We then invoke results from our previous work on the 
hamiltonian formulation of the FQHE [|l~0| p~2[| which al- 
lows one to rewrite the hamiltonian in terms of composite 
fermion degrees of freedom. The advantage is that these 
particles see an effective field which is just right to fill 
an integer number of LL's, thus yielding a nondegener- 
ate ground state on which we can base a Hartree-Fock 
calculation of gaps. 



Specifically, the hamiltonian approach provides an ex- 
pression for the LLL projected electron density p(q) in 
the CF basis JToj— . Since the effective hamiltonian is 
expressed in terms p(q), this all we need. 

In Section II, we derive the effective hamiltonian and 
in Section III we express it in the CF basis. Section IV 
contains the calculation of gaps and conclusions follow in 
Section V. Details are relegated to the two Appendices. 
While this paper focusses on the calculation of transport 
gaps, the approach presented here can be extended to cal- 
culations of other physical quantities in the Hamiltonian 
theory ||. 



II. THE EFFECTIVE 
HAMILTONIAN-ELECTRONIC BASIS 

Our starting point will be the full electronic Hamilto- 
nian, which is given by 

H = ^2 ala a Ei + - / ^pi 2 (q)p34(-q)w(g)a t 1 4a 4 a 2 (1) 

= H a + V (2) 

where any subscripts a stands for pair of indices 
(n a , m a ) (n labels the LL index and m the angular mo- 
mentum within the LL) . Both labels assume integral val- 
ues from zero to infinity and 



E a = n a oJo 



(3) 



with loq being the cyclotron frequency. We have 
dropped the zero-point energy and kept only the normal- 
ordered part of the interaction. The operators a and 

are electron destruction and creation operators, v(q) 
is the electron-electron interaction and p a p{<\) are the 
matrix elements of the charge density operator. These 
matrix elements can be determined using the standard 
decomposition of the electron coordinates and momenta 
into cyclotron (rj e ) and guiding center (R e ) variables 
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(4) 
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where II e = p e + eA is the velocity operator of the 
electrons and I = y/l/eB is the magnetic length. These 
coordinates obey the commutation relations 



\jlex -i Ve 



-iF 



[Ve R e 



■■ 
(6) 



Expressing the electron coordinate as r e = R e + r] e we 
have 



Pa/3(q) 



-iq-(R c +TJ e 



)\n p mp) 



m Q e 



-iq-R E 



(m) (n) 

Pj x Pj 



(7) 

Imp) x {n a \e-^'h\np) (8) 

(9) 



where the superscripts m and n indicate that the ma- 
trix elements correspond to the guiding center and cy- 
clotron coordinates respectively. These matrix elements 
can be expressed in terms of Laguerre polynomials. How- 
ever, for this section we only need to know that 



(™) i \ 

Poo (q) = e 



/4 (10) 

^P$\*)p%\*) = P%\*+*)e^*« fi ' 2 (11) 

& 

which reflects the magnetic translation algebra: 



Here we have removed the normal-ordering, which 
leads to an unimportant additive constant in the energy. 
The operator p(q) has the same commutators as e 4q R . 

To extract the leading correction due to higher LL's 
we write the Schrodinger equation as 



Hoo H 0n > 







= E 


~ 4>~ 











(17) 



where <j) is restricted to the space spanned by Fock 
states composed of just the LLL states and \ stands for 
everything else. The exact equation obeyed by 4> is 



Hoo + H[ 
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E — H n i n i 



H n , )</> = E<t> 



(18) 



which is not an eigenvalue problem since E appears on 
both sides. However we may approximate as follows: 
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E - H n ' n ' 



Hn'n> 



+ 0{v/w) 



(19) 



since the eigenvalue E we are interested in is of order 
v and the eigenvalues of H n ' n i arc of order luq. To the 
same accuracy in k = v/loq we can also replace 



H n ' n ' — H n , n , - ^2 a) a a a n a 



UJo 



(20) 
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It will prove convenient to define 



f f d 2 q 



(12) 



(13) 



where the potential is included as part of the measure. 

The leading term in the effective hamiltonian is the 
obvious contribution from Hqo, the restriction of H to 
the LLL: 

Hoo = E fl Om fl 0m ■ ■ U>0 
m 

1 f r t t 

+ 2 J y a 0m 1 a 0m 3 a 0m i a0m 2 

= ^<«l3 a -4a m2 ^ ) m2 (q)^ 3 U(-q) e " 92 ' 2/2 (I 4 ) 

where we have labeled LLL operators with just an m 
index. In terms of p(q), the density restricted to the LLL 
and defined by 



In fact, one can do a little better than this, and in- 
corporate the interactions into H n r n i in by including the 
Hartree-Fock energy of the n' th level. This can be viewed 
as a perturbation theory in the non-HF part of the in- 
teraction. We will present results both with and without 
this HF energy in H n i n i. To maintain generality we de- 
fine a reduced energy 

e(n)=n+^-J ^ v{q )(\ p $ (fif - (21) 

where we have used the fact the the FQH states have 
uniform occupations of v in the LLL, and we have nor- 
malized the LLL energy to be zero. This leads us to the 
effective Hamiltonian 



n 00 



H 00 + SH t 



00 



where 
6H = 



H 



On' 



H n '0 



a l a 6 a 5 a 7 
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(22) 
(23) 
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H°, 

n'n' 



■a\a) 3 a 4 a 2 



P(q) = ^2aL i a mj Pm t m J (q) 



we can write 
Hoo = 



V= 2/ e 



-{qlf/2 



p(q)p(-q)- 



(15) 



(16) 



'<?l J 92 

Pl2(qi)P34(-qi)P87(q2)P65(-q2)] (24) 



The prime on Y^' denotes the fact that if SHqo is to 
act on the LLL sector, we need 

n$ = n e = 114 = n 2 = 0. (25) 

This means we can have two kinds of terms: 
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• 5Hqq , a term in which both a\ and a\ excite the 
electron to higher LL and 07 and 05 bring it back 
to the LLL. 



G(qi , q 2 ) p™ (-qi)/4"l 4 (-qi)pS } (-qa)p^U (-<fe) 



a m s a m 6 a m 5 Ct m3 a m 4 a m 2 ] 



(31) 



• SH^ in which either one of a\ or ag excites the 
electron to higher LL and one of 0,7 or a.5 brings it 
back to the LLL. 



,t 



It is readily seen that 



-i( qi Xq 2 )! 2 /2 -(qr 2 + g 2 )Z 2 /4, 



(qi + q2)Pm3m 4 (-qi)/0m 6 m 5 (-q2) 



a m 8 a mn a m 5 a m 3 a m 4 a m 2 ] 



m(2) _ 

dii 00 - 



91 -'92 



agaga4a 2 - 



X!'° 12 (q i )^34(-qi)p87(q2)/365(-q2) G(qi,q 2 )= 
1 



1... 



ri'>0 



Po» ) >(q2)p.^o(qi) 

e{n')uj 



(32) 
(33) 



-(^53^71 - ^51^73) 



' (e(n[) + e(n' 3 ))uJo 
The contribution from the first delta function is 



(26) If we now rearrange the operators (whose subscripts 
now all describe m) so as to match those on p.y as follows 
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(34) 



c,/ x ^ Po'l(q2)^o(qi)po?n 3 (-q 2 )Pn3 ; ,o(-qi) , 07 ^ r m 

^il9ij92j— 2_^i 7 — TT — 7 — \ \ 'Jwe arrive at the following expression for SHq ' 
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ni ,713=1 



e(m) +£(^3) 



Using completeness to do the sum over mi and 7713, < 
using the magnetic algebra Eqn. (|ll| ) and doubling the 
answer because of the other delta function which makes 
an equal contribution (upon relabeling of dummy indices) 
we arrive at 



(i) 
on 



G(qi,q 2 )e 

qi -Jqi 

i(qiXq 2 )( 2 / 2 



-( 9 2 + 9 2 )Z 2 /4 



p(qi +q 2 )p(-q 2 )p(-qi) 

P(q 2 )p(-q 2 ) - e- l(qiXq2)i2 p(qi)p(-qi)" 



(35) 



SH, 



(2) 
00 



-1 

2uj 
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<?l<32 



-iqixq 2 ; 2 ^/„ ,„ \=/ . _ \ ip / „ \ In writing the answer m terms of p s we must pay 

p(qi +q 2 )p(-qi-q 2 )^i(<7i,9 2 ) ,. , °, , . , K , , 

attention to tne order since tney do not commute with 
(28) eacn other. It should also be noted that the first term 
of Eqn. (H|) produces a contribution to another effective 
two-body term via 



Now for the term with just one higher LL excitation: 



5H {1) - -- 



,t„t. 



[P12 (qi V34 (-qi ) ps7 (q 2 ) pes ( - q 2 ) 



p(qi + q 2 ) 



27T/ 2 



(27r) 2 (5 2 (qi + q 2 )+ : p(qi + q 2 ) 



,t„t. 



(29) 



Consider the contribution when a\ takes the electron 
up to a higher LL and and 07 brings it down. The other 
three possibilities obtained by permuting 1 < — ► 3 and 
5 < — > 7 make the same contribution upon relabeling of 
dummy indices. The net contribution due to these terms 
in which a single excitation to a higher LL occurs is 



where the first term contains the density v/2~kI . 
The final effective hamiltonian is given by 



H, 



eff 



00 



SH 



(2) 



SH, 



(i) 



(36) 



(37) 



where Hqq, SHqq and SHqq are given by Eqns.(fL4| 
|) and @. 



SH$=- / Ya\n B a L & arn 5 a m7 S7ia\ ni a\ n3 a mi a m2 

J ai J m. 1 o 

/034(-qi)p6s(-q 2 ) (30) 



qi Jqi 1. 

psi(q 2 )pi 2 (qi) 



e(ni)u! 



Pm 8 m 1 (q 2 )Pm 1 m 2 (qi)x 
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III. THE EFFECTIVE HAMILTONIAN IN THE 
CF BASIS 

The preceding calculation of the effective LLL theory 
in terms of electrons could have been performed decades 
ago since it is not dependent on any particular repre- 
sentation of p. However, written in terms of electronic 
variables, the effective theory suffers from the same hand- 
icap as the original one: the noninteracting part of the 
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hamiltonian has a huge ground state degeneracy for fill- 
ing fraction v = 2p s+i < 1- To get a nondegenerate 
starting point, we must transcribe H e ff(p(q)), calculated 
above in the electronic basis, to the CF basis. The CF 
description will be in terms of particles that see a weaker 
field of just the right strength that they fill exactly p 
LL's, with no degeneracy in the ground state. Over the 
years we have developed a route that takes one from 
the electron-representation to the CF-representation. We 
will not describe it here since it is long and has been de- 
scribed before ]To|-[l^|. All we need here are the end 
results, which look a lot more attractive than the path 
that led to them, fraught as it was with approximations 
and inspired guesses. Instead, we will arrive at our final 
result armed with hindsight. 

Consider a CF Hilbert space where each fermion is de- 
scribed by a coordinate r and momentum p. (Hereafter, 
all unsubscripted variables will represent the CF degrees 
of freedom.) From these we construct the velocity oper- 
ator 



n 



eA* 



.4* 



2ps + 1 



(38) 



where the weakened vector potential A* is what the 
CF sees. In terms of these variables, the electron guiding 
center R e takes the form 



R r 



I 2 



(1 + c) 



-z x IL 



2 2 PS 



2ps + 1 



(39) 
(40) 
(41) 



Note that R e obeys the right commutation relations. 
Also, since it is written in terms of an object that sees 
a weaker field A* — A/(2ps + 1), if we express H e ^ 
in terms of these variables, we will not encounter the 
degeneracy problem at the Jain fractions. This is the 
motivation for switching to the new coordinates. 

It is useful to write R e in terms of the CF guiding 
center and cyclotron coordinates R and 77: 



R, 



R + crj. 



(42) 



Since we have embedded R e in a regular fermionic 
Hilbert space, we have room for another pair of guiding 
center- like coordinates. Let us call them R„. These will 
be of the form 



R„ = aR + (3r\ 



(43) 



Demanding that R„ commutes with R e gives us — 
a/c. How about the overall scale of the operator? Here 
is where our previous work tells us to choose a = 1, that 
is 



R„ = R + 77/c. 
In terms of r and II 



(44) 



R„ = r 



I 2 



c(l + c) 

[Rvx ) Rvy\ — ^ / C , 

[R e , Rt,] = 0. 



z x n 



(45) 

(46) 
(47) 



The merit of this choice is the following. These com- 
mutation relations correspond to the guiding center co- 
ordinates of a particle of charge — c 2 = —2ps/(2ps +1). 
This is precisely the charge of an object that must pair 
with the electron to form the CF and we refer to it as the 
pseudovortex coordinate, since it has the same charge as a 
2s- fold vortex in Laughlin states. It must be emphasized 
that R v cannot be directly identified with the physical 
vortex that exists around an electron. There is, how- 
ever, a connection between R„ and the physical vortex 
which arises upon the choice of a suitable state in the R„ 
variables, as explained below. 

R u is a cyclic coordinate that does not enter the hamil- 
tonian, and so has no dynamics. The eigenfunctions of 
H will be of the form 



^(R ex ,R vx ) = V{R ex )^(R vx ) 



(48) 



where we have chosen as a commuting pair of coor- 
dinates R ex and R vx with conjugate momenta R ey and 
R vy . (One can also use the Bargman representation and 
consider ^(z^), where z and w are complex numbers, 
as we will find convenient to do later.) Whereas ^(R ex ) 
is an eigenfunction of H(R e ), ^(R^x) is completely ar- 
bitrary. Thus each eigenfunction is infinitely degenerate. 
This is exactly the kind of degeneracy a gauge symme- 
try would introduce. We must therefore "fix our gauge" 
i.e., choose an arbitrary function ^>q(R vx ) to accompany 
*i>(R ex ). No physical observable (function of R e ) will 
depend on this choice. Our previous derivation based 
on canonical transformations [p~0| — naturally led to the 
following constraint on ^>(R VX ) 



^2 e -iq Rvj j (physical states) = 



(49) 



The constraint arose because we had introduced addi- 
tional oscillators at the cyclotron scale and these had to 
be paid for. 

In hindsight we see that we are free to take an ax- 
iomatic view of R e and R v and to simply introduce 
them as the above functions of CF coordinates and mo- 
menta (Eqs. ( p9[^5| )). The hamiltonian depends on R e 
alone. The constraint we were led to in our earlier work, 
Eqn. (|4S|), is an acceptable choice for fixing the gauge. 
However, since all physical observables depend on just 
R e , and R e commutes with R„, we may assign to R„ 
any dynamics we want without changing the physics. 
For example, we could add to i?(R e ) a piece H(R, V ) 
which is any generic repulsive interaction and demand 
that ^o(R vx ) be its ground state. Now it turns out that 
for any generic repulsion there is essentially a unique an- 
swer, the v = l/2s bosonic Laughlin wavefunction p4[. 
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The reason is as follows. First, we must choose R v to be 
a bosonic coordinate (the wavefunctions have to be an- 
tisymmetric in r e and R e , and hence symmetric in R„). 
Next, since has a magnetic algebra charge of — c 2 , and 
there is one pseudovortex per electron, it is easily seen 
that it is always at filling 



P 



2ps + 1 



2ps + 1 2ps 



1 

2s' 



which leads us to the wavefunction 



Zi) 2s e 



(50) 



(51) 



Note that the magnetic length appropriate to the bosons 
(with charge c 2 e) has been used. 

Now it might seem that all these virtues are moot since 
H v is a fictitious coordinate whose dynamics should make 
no difference to the actual physics. While this is certainly 
the case in any exact calculation, in the approximate 
HF calculations we will employ (where the wavefunction 
does not have the factorized form of Eqn. (f48| ) , different 
choices give different answers and the above choice gives 
the wave functions with the best correlations. 

To see this, let us first define the CFHF wavefunc- 
tion. Here we work in the CF single-particle basis in 
the effective field B* (defined in more detail in the next 
section), and fill the lowest single-particle states upto a 
filling of v* — vj{2sv — 1), which reduces to v* = p for 
v = pj (2ps + 1). It is easy to show |ll| that these states 
are indeed HF states of the Hamiltonian. 

As mentioned above, the CFHF wavefunction cannot 
be a true eigenstate, since it not of the required product 
form of Eqn. (|4§|). However, we can "project" it in the 
following way 

Vy C FHF({~Re, l ,~RvA) = (52) 

|*„({R^}) >< *„({R„ !l })|* CFffF ({R e , i ,R lVi }) > (53) 

= *«({R«, < })*e,*.({R e ,i}) (54) 

The result is a wavefunction which is explicitly in the 
form of a product. Clearly, different choices of ^ v will 
lead to different results for the electronic part of the wave- 
function, and this dependence is indicated by the sub- 
script on the electronic wavefunction. For a given CFHF 
wavefunction, say that of the ground state, we can imag- 
ine systematically varying to obtain the electronic 
wavefunction with the best possible energy. It would 
then be reasonable to call this ^f v the best possible pseu- 
dovortex wavefunction. 

In Appendix A we present the details of how to carry 
out the "projection" . The result can be gleaned from an 
identity which we borrow from Appendix A; 



d 2 z v f({z Vli })e 



£■ 



-£■ 



/(KJ) 

(55) 



which shows that correlations of the wavefunction 
get transferred to the ^ e wavefunction upon "projec- 
tion" . Thus, the zeroes of the bosonic wavefunction for 
the pseudovortices induce Laughlin- Jastrow factors in the 
electronic wavefunction. This establishes the precise con- 
nection between R„ and the physical vortices. 

Furthermore, it is well-known that wavefunctions with 
2s Laughlin- Jastrow zeroes attached to the electrons 
have very nearly the best energy for generic repulsive 
electronic interactions at v = p/(2sp + 1). Thus, the 
ty v which produces these correlations upon "projection" , 
namely the v — bosonic Laughlin wavefunction (Eqn. 
(pl|)), must therefore be very nearly the best choice for 



IV. DETAILS OF THE HARTREE-FOCK 
CALCULATION 



Turning to the HF calculation of gaps we will once 
again use the preferred density []l3|| in H as we have in 
the past: 



p(q) - c 2 x(q) 



(56) 



This should make no difference in any exact calcula- 
tion of gaps, which can be seen as follows: First note 
that when we expand out H p = H(ff), there will be a 
p(q)p(q) term, which is there to begin with, a x(q)x(q) 
term which contributes a constant depending on our 
choice of (equivalent to a choice of "gauge"), a con- 
stant which will drop out in the difference between the 
ground state and a state with a widely separated particle- 
hole pair, and finally a p(q)x(q) part whose expectation 
value will vanish as along ty v is translationally invariant. 
The reason we make the replacement p(q) — ► p(q) p is as 
in our earlier papers: in the HF calculation, this choice 
builds in Kohn's theorem JlH] (leading to S(q) ~ q 4 fli6| ) 
as well as the correct charge and dipole moment Jl7| of 
the CF at tree level, thus making it plausible that the 
results do not suffer strong vertex corrections pj| . 

It is now straightforward, though tedious, to compute 
the correction to the transport gap from the terms in- 
duced in the effective Hamiltonian. One follows the rules 
of standard first order perturbation theory and takes the 
average of the perturbation in the unperturbed states. 
In this context, note that the HF states of Hqo are not 
necessarily HF states of Hqo '■ However, this effect is 
of order v(q)/u>o, and will affect the gaps to second or- 
der in v(q)/u)Q. In order to evaluate the averages in the 
HF states, one needs to do two two-dimensional integra- 
tions. While it is possible to do these numerically for an 
arbitrary potential, we have simplified the problem by 
choosing a particularly tractable form for the potential 



v(q) 



2ve 2 



2ve 2 



(57) 
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This potential enables us to perform all the integrations 
analytically. The relevant integrals are tabulated in Ap- 
pendix B. 

Let us briefly sketch how the process works. We first 
represent the projected preferred density in terms of CF 
annihilation and creation operators (d, d') in the Landau 
gauge: 



ipCF(r) = E 4>n,xd n ,x 



(58) 



n.X 



4> n<x = cgi^/cnvc*- v 2 / 2 ^ 2 H n (( x - x)/i*) (59) 



where C is a normalization constant, I* — l^/2ps + 1 is 
the magnetic length in the effective field, and the H n are 
the Hermite polynomials. In terms of d ni x and d' n x , the 
preferred density looks like 



p(q) = £ e^ x d} Qyl *d Qyl *p ni7 Jq) (60) 

£ — * ni,X % — ™2,^M 2 — 

{n,}X 



where the matrix element is defined by 

P ni nM = (-l)»<+"3y^e^-^)(<W 2 ) 
\cQ/V2) lni - n2l e- c2 Q 2 / 4 Ltr n2l (c 2 Q 2 /2) 
~ c 2 (Q/cV2)^- n ^e-Q 2 ^ c2 L l ^- n2l (Q 2 /2c 2 ) ) (61) 



The two-body and three-body terms have to be treated 
differently. The two-body terms all involve the matrix 
elements of p(q)p(— q) integrated against different mea- 
sures. One first computes the contribution to the gap due 
to /o(q)p(— q). Recall that the transport gap for v = | is 
the energy difference between the n = CF-LL and the 
n = 1 CF-LL. Thus, 



S(q) = <n= l\p(q)p{-q)\n = 1 > - 
< n = 0\p(q)p(-q)\n = > 

This quantity is easily found to be 



(62) 



l-c 2 



-c 2 (ql) 2 /(2(l-c 2 )) 



- 4 C 2 + _2£^(^)2j e -(c 2 +c-^)( 9 2 /(4(l-c 2 )) 

+ ^2C 4 - T £l 7 (g/)2^ e -(«0 2 /(2c 2 (l-c 2 )) 
+ ( (? /)2( 1 _ c 2) e -(l-c 2 )M) 2 /(4 C 2 ) 



in the n = and the n = 1 CF-LL's. After some effort, 
the expectation value of this operator in the state |m > 
can be computed to be 

< m\ : p(qi + q 2 ) : p(-qi)p(-q 2 )\m > = 



C'2 



qixq 2 (r) 2 



E (l-iV F (m))(l-iV F (n 2 )) x 



(-qi)p it 2 in 

(-02 j 

. e -|qixq 2 (i*) 2 N F (ni))N F (n 2 ) X 

Pm ni (qi + q2)p ni „ 2 (-q2)p„ 2m (-qi) 
e-jqixqaC*) 2 £ AT F (m)(l - N F {n 2 )) x 

PmmC-qOPm^Cqi + q2)Pn 2 m(-q2) 

_ e § qi xq2(r) 2 £ (i_ jv F (ni))iV>(n 2 ) x 

ni"2 

Pmn 1 (-qi)Pn 1 ™ 2 (-q2)Pn2m(qi +qa) 

+ e -iqixq 2 (r) 2 N F (ni)N F (nQ) x 

rtin 2 

Pmn 1 (-q2)p„ 1 „ 2 (-qi)^n2m(qi + qa) 



qiXq 2 (r 



)2 E ^.(nOCl-JVpfna)) x 



Pmm (-q2)P ni r>2 (Qi + q2)Pn 3 m(-qi) 



(65) 



Only the n = CF-LL is occupied, and the infinite 
sums over m and n 2 (where not truncated by N F ) can be 
carried out using the following expression for the matrix 
elements 



P, 



lin M =< mle-^ - c 2 e~^' c \n 2 > (66) 



and using the completeness of the set of states \n>. Suf- 
fice it to say that the results of this calculation can also 
be expressed entirely in terms of powers and gaussians. 
Appendix B describes an integral that can be used to 
find both the two-body and three-body contributions to 
the gap. 



V. RESULTS 



(63) 



It is seen that the result contains only powers of ql 
and gaussian factors. With the assumed expression for 
the interaction, Eqn. (57), the integrations can be carried 
out analytically. 

The contribution of the three-body term to the gap is 
somewhat more involved. It can be seen that it involves 
the expectation value of the operator 



: P(qi + q2) : /9(-qi)p(-q2) 



We have carried out this calculation for three cases, | 
and | (which are spin-polarized), and spin-singlet |. The 
common factor in all these cases is that only the n = 
CF-LL is occupied. The only difference is in the value 
of c, which is c = \/2/3 for | and \/4/5 for both | and 
singlet | . This leads us to the interesting fact in our CF- 
HF approximation there is no difference between ^ and 
singlet |! In computing the effect of LL- mixing excita- 
tions to all electronic LL's must be included in principle. 
In practice, we found that the contributions drop rapidly 
for high n'. We found that a cutoff of n' max = 30 was 
sufficient to capture all the contributions to five-figure 
accuracy. 

In Figure 1 we present the variation of the spin- 
polarized and spin-reversed transport gaps for | as a 
(64) function of k — (e 2 /el)/oj c for a "thickness parameter" 
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0.12 



0.1 - 



0.08 - 



A 0.06 



0.02 













A with HF energy 




A SR with HF energy 




A without HF energy 




A SH without HF energy 



0.03 



0.02 - 



0.01 





*•• • . 




A with HF energy 




A SR with HF energy 




A without HF energy 




A SH without HF energy 



FIG. 1. Variation of the spin-polarized and spin-reversed 
gaps with LL-mixing for v = i at A = \.%. The solid and 
dotted lines have the electronic HF energy included (Eqn. 
(Ell)) while the other two do not. All energies are in units of 



0.12 




FIG. 2. Variation of the spin-polarized and spin-reversed 
gaps with A for v = i at k = 2. For reference, the results 
without LL-mixing are also presented. All energies are in 
units of e 2 /el. 



A = 1.21 (for which the HF gap in the limit of no LL- 
mixing has appromimately the same value as for the true 
gap for Coulomb interaction). The electronic HF energy 
(Eqn. (^l|)) has not been added to the cyclotron energy 
in the dashed curve (which therefore shows a linear de- 
pendence on k), while it has been included in the solid 
curve. It is seen the the gap decreases by a few percent 
for realistic k, but there is no dramatic effect. 

In Figure 2 we show the variation of the gaps with A 
for fixed k = 1, this time also plotting the result in the 
absence of LL-mixing. It is seen that high values of A 
suppress Landau level mixing. These results are in agree- 
ment with the previous fixed-phase diffusion Monte Carlo 
(FPDMC) studies || at a semi-quantitative level. In gen- 
eral, our computation shows the gap to be more robust 
than the FPDMC calculation. For instance, at k = 4. 



FIG. 3. Variation of the spin-polarized and spin-reversed 
gaps with LL-mixing for v — | and spin-singlet | at A = 0.8/. 
Note that the spin-reversed gap applies only to | while the 
spin-polarized gap applies to both | and singlet t. The solid 
and dotted lines have the electronic HF energy included (Eqn. 
(Ell)) while the other two do not. All energies are in units of 



which roughly corresponds to r s = 10, the FPDMC re- 
sults for the pure Coulomb interaction show a reduction 
of about 30% in the spin-polarized gap || for | . In con- 
trast, our results for A = 1.2, which has roughly the same 
gap in the absence of LL-mixing, shows a reduction of 
only about 15% (with the electronic HF energy included 
in H n / n i). On the other hand, this reduction is about 
the same as the FPDMC results for a Coulomb interac- 
tion corrected for sample thickness (thickness parameter 
= 1.51) §. Our calculations for A = 1.21 clearly do 
include the strong suppression of the interaction at large 
wavevectors, which is characteristic of the effects of sam- 
ple thickness. However, in order to carry out a detailed 
quantitative comparison one would have to re-do our cal- 
culation with the same interaction potential as was used 
in the FPDMC work. In the present work, we have con- 
tented ourselves with a proof of principle, concentrating 
on the analytically tractable potential of Eqn. (|57|). 

Figure 3 displays the variation of the gaps with k for 
■i/singlet | (note that the spin-reversed gap is only for 
I), this time for A = 1.6Z (which again gives a LLL gap 
roughly equal to that obtained from exact diagonaliza- 
tion for the Coulomb interaction for v = g). 

Qualitatively, it is seen that the FQH states are very 
robust to LL-mixing, at least as far as the transport gaps 
are concerned (which has been pointed out before ||,|8|Q ) ■ 
While the calculations which leave out the electronic HF 
energy in H n > n ' show a linear dependence on k, the ones 
which include the HF energy show a much more physical 
behavior, with results that are even more robust under 
LL-mixing. 
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VI. CONCLUSIONS 

The FQHE has been a fertile source of new ideas in 
the last two decades. Much of the initial insight into the 
physics of the FQHE was obtained by considering wave- 
functions ||^|]. However, calculations of physical quanti- 
ties based on the wavefunction approach were forced to 
resort to computationally demanding techniques. The 
advent of Chcrn- Simons field theories §J-f| of the 
FQHE raised the possibility of approximate analytical 
schemes to calculate the desired properties. Our Hamil- 
tonian approach lies squarely in this class of approaches. 
Building on previous results, we were able to obtain a 
LLL theory of CF's. Simple approximations such as HF 
in this theory produce results [ p"3| in very good agree- 
ment with those obtained numerically on small systems, 
and with experiments in high fields. However, as the field 
is lowered, the question of LL-mixing becomes unavoid- 
able, and must be accounted for theoretically. Previous 
approaches to this issue have been based on exact di- 
agonalizations Q or the Quantum Monte Carlo method 
on finite systems. 

The approach to LL-mixing presented in this paper 
is analytical and approximate. Despite the approxima- 
tions, the results agree at a semi-quantitative level with 
previous results |(],|8|||. For illustrative purposes, we have 
chosen to compute the effect of LL-mixing on the trans- 
port gap. In principle, one can calculate the effect of LL- 
mixing on any physical property in this approach. Some, 
such as dynamical response functions, or finite tempera- 
ture spin polarization, would be prohibitive to calculate 
numerically, but are easily computed in our approach 
fl3|| . One interesting question concerns instabilities of 
the FQH states to LL-mixing. This could be addressed 
in our approach by applying a conserving approximation 
to calculate the collective modes ^,^^3), and seeing 
whether any instabilities develop in these modes as LL- 
mixing increases. We intend to pursue this and other 
questions in the future. 

We are grateful to the NSF for partial support under 
grants DMR0071611 (GM) and DMR0103639 (RS). 



VII. APPENDIX A 

In this appendix we will examine how "projection" of 
a CFHF wavefunction to a product form works in detail 
for v = |. In this case the ground state CFHF wavefunc- 
tion is particularly simple, being just the wavefunction of 
a fully filled CF-LL. In order to carry out the projection 
we first have to reexpress this wavefunction in terms of 
{R e ,i, Rt, : i}. Let us focus on the single-CF states in the 
lowest CF-LL. These states have the property that they 
are annihilated by the CF-cyclotron destruction opera- 
tor. The CF-cyclotron operators can be expressed as 



R, 



R 



CF,x 



iR 



CF,y 



2 



(67) 



CF,x 



R 9 



iR C CF,y — I 

^CF,y = ' 

iR 9 — I 

lJX CF,y ~ 1 



2 *ot 



1-c 



; 2 *> at 
l— c 2 g 



where 



[a c ,al] = [a g ,a\} = 1 
[a g ,al] = 



\a c , a 



tl - 



(68) 
(69) 
(70) 

(71) 
(72) 



We can now express these operators in terms of electron 
and pseudovortex guiding center creation and annihila- 
tion operators, which are defined by 



R e>x + iR e , y = lV2A\ = f - 2l 2 8 e 

Re,x ~ iRe,y = lV^A e = f + 2l 2 d e 
hH A — hL _l 2fn 

c S± v — 2 -I- c 2 U v 

lV2 4f _ £^ _ 2l 2 a 



Rv : x 

R'D ,X 



iR 



v,y 



(73) 
(74) 
(75) 
(76) 



where d e , d e are shorthand for etc. The result- 

ing expressions are 



(77) 
(78) 
(79) 
(80) 



,t _ 





cAl 




-c 2 


At- 


cA v 


Vi 


-c 2 


A v - 


cAl 


Vi- 


-c 2 


Al- 


■cA a 


Vi- 


-c 2 



Consider now the CF state with CF-LL index and an- 
gular momentum zero. This state will be labelled |0,0 > 
as satisfies 



a c |0,0 >= (Ay - &4)|o,0 >= 
a 9 |0,0>=(A e -c4)|0,0>=0 



(81) 
(82) 



Since all the available states are tensor products of LLL 
states of electrons and pseudovortices we can express this 
state as 



k,,k v =0 



|0,0>= A(k e ,k v )\k e > e ®\k v > v (83) 

where 



I, I 



■|0> e 



y , 27r2 fc <!fe c ! V ' 



-\*.\V*F 



I ky 



(Al) k v 



|0>„= 



p - C 2 | Z „| 2 /4i 2 



(84) 



(85) 



are the LLL states of electrons and pseudovortices. 

Applying the conditions of Eqs( |8l| , p2|) it can easily be 
found that in 10, >, one is forced to have k e = k v and 



|0,0>= Vl^? ^ c k \k> e ®\k > v 



2-kI 2 



(86) 
(87) 
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By similarly using the expression for in terms of the 
electron and pseudovortex coordinates and their deriva- 
tives one can verify that 



(at) m 
|0,m>= V 9> 



|0,0>= 



1 



T) |0 < 0> 

(88) 



Thus, filling all the CF states in the CF-LLL leads to the 
wave function 



, op=1 = CU^e,i-Z e , j )\0,0> 



i<j 



z e ,j)e~ 



- 



72— + " 



(89) 
(90) 



where C is a normalization constant. It can be seen that 
this is not in the form of a product of functions of z e and 
z v alone. 

Now consider "projecting" this wavefunction into a 
product form against an arbitrary wavefunction *&„ de- 
pendent on {z v ^}. Performing a few gaussian integra- 
tions one can easily show that 



„ _ — ip — 

d z z v f({z v ,t})e " 



(91) 



Now if one takes to be the v — | bosonic Laugh- 
lin wavefunction, it can be seen from eq(|9l]) that the ze- 
roes of the pseudovortex wavefunction are translated into 
Laughlin-Jastrow factors in the electronic wavefunction. 
This establishes the precise connection between correla- 
tions in the pseudovortex coordinates and correlations in 
the projected electronic wavefunction. 

VIII. APPENDIX B 

In this appendix we present an integral that can be 
used to find the contributions of both the two-body and 
three-body terms. It is found that the contributions to 
the gap in both cases can be expressed as a product of 
powers of (qil) 2 , (qzl) 2 , or (qi + q 2 ) 2 l 2 , and gaussians 
of the same arguments. This has to be integrated with 
a complicated measure in qi and q2. Calling q\l = y\ 
and q 2 l = 2/2, and using the specially chosen form of the 
interaction, it is easy to see that all the integrals can be 
reduced to special cases of the following integral 

2-k oo 

Qn({3,a 1 ,a 2 ,a 3 ,m 1 ,m 2 ,m 3 ) = J 4§^r / dy 1 dy 2 

o o 

izlT-^ffiJ^^ y™ly™2 e im 3 (02-0l) e ^ 2 (l+2/3)qiXq 2 (92) 

One can perform the angular integration by choosing 
6 = 82 — #i, and using the identity 



(16 
2^ 



(t \ n 
-J / W (2V56) (93) 



where /„ are the Bessel functions with imaginary argu- 
ment. 

This identity holds even for negative and/or complex 
values of a and 6, which for our integral attain the values 
a=3 + f-f and6=-i-|-^. 

For the radial integrals one uses the following two for- 
mulas J27J (in which J v are the Bessel functions, $ is the 
confluent hypergeometric function, and F is the hyper- 
geometric function) 



00 2 

J x^e- ax J v {fix) 





2"+!q — 1 — r(v+i) 

-stib—l; 



*( V + l i + 1 ;v + l;-p 2 /4a) (94) 



/ dte- st t b - 1 ^(a; c; kt) = r(b) S - b F(a, b; c; k/s) 

[\s\ > \k\] 

= T(b)(s - k)- b F(c - a, b; c; k/(k - s)) 

[\s-k\ > \k\] 

to obtain, for the case ab > 
Q n {(3,a 1 ,a 2 ,a 3 ,m 1 ,m 2 ,m 3 ) = L^l 2 - 1+ ^+ m ^/ 2 x 

|-OMm.s+n r(n+(mi+Tn3 + l)/2)r(7i+(m 2 +m 3 + l)/2) x 

f(Ti+(mi+m3 + l)/2,n+(m 2 +m 3 + l)/2;n+m 3 +l;z) 
(l+2A 2 +2a 1 )™ + (™i+™3 + i)/2(i +2 A 2 +a 2 )" + ( m 2 + '"3 + i)/2 



(95) 



(96) 



where 



Aab 



' (1 + 2A 2 + 2ai)(l + 2A 2 + a 2 ) 
Similarly, for the case ab < we obtain 



(97) 



Q n (P,a 1 ,a 2 ,a 3 ,m 1 ,m 2 ,m 3 ) = L^l 2 m ^ 1 (2b) m3+n x 
(-4a&)- rl ^ m2+m3+1 )/ 2 (i + A 2 + ai )^- m ^/ 2 x 

r(n+(mi+m 3 + l)/2)r(n+(m2+m 3 +l)/2) n+ ( TO2+m3+ l) / 2 x 
n!()i4-m 3 )! 

F ( m 3 -m 1 +1 ^ n + m 2 +m 3 + l ;n + T n 3 + l;z) 



(98) 



where 



-Aab 



-4ab + (1 + 2A 2 + 2ai)(l + 2A 2 + a 2 ) 



(99) 
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